Magnetic penetration-depth measurements of a suppressed superfiuid density of 
superconducting Ca .5Na sFe 2 As 2 single crystals by proton irradiation 
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We report on the dramatic effect of random point defects, produced by proton irradiation, on 
the superfiuid density p s in superconducting Cao.5Nao.sFe2As2 single crystals. The magnitude of 
the suppression is inferred from measurements of the temperature-dependent magnetic penetra- 
tion depth A(T) using magnetic force microscopy. Our findings indicate that a radiation dose of 
2xl0 16 cm -2 produced by 3 MeV protons results in a reduction of the superconducting critical tem- 
perature Tc by approximately 10%. In contrast, p 3 (0) is suppressed by approximately 60%. This 
break-down of the Abrikosov-Gorkov theory may be explained by the so-called "Swiss cheese model" , 
which accounts for the spatial suppression of the order parameter near point defects similar to holes 
in Swiss cheese. Both the slope of the upper critical field and the penetration depth A(T/T c )/A(0) 
exhibit similar temperature dependences before and after irradiation. This may be due to a combi- 
nation of the highly disordered nature of Cao.5Nao.5Fe2 As2 with large intraband and simultaneous 
interband scattering as well as the s ± -wave nature of short coherence length superconductivity. 



I. INTRODUCTION 

Proximity of the superconducting and magnetic states 
in iron-based superconductors has stimulated exten- 
sive studies of the gap nature^— order-parameter 
symmetry^— and the pairing mechanisms in these 
materials.— The response of the superconducting conden- 
sate to impurities is sensitive to the symmetry of the su- 
perconducting state, and their influence has been widely 
investigated to gain better understanding of the nature of 
the order parameter in both low- and high-temperature 
unconventional superconductors^— 

The Abrikosov-Gor'kov (AG) theory^ 3 - explains the ef- 
fects of impurities in the low-T c superconductors, where 
a large superconducting coherence length £ effectively av- 
erages the suppression of order parameter at the impurity 
sites over many impurities, leading to a uniformly sup- 
pressed order parameter. However, the AG theory breaks 
down when applied to the effect of disorder on super- 
conducting properties in the cupratcs superconductors^ 
where £ is short and comparable to the average spac- 
ing between disorder centers. The order parameter is 
therefore suppressed locally at the impurity site and has 
a chance to recover between impurities. The influence 
of disorder on the superfiuid density p s in cuprates is 
well described by the so called "Swiss cheese" model, 
which considers spatial dependence of the order pa- 
rameter and its strong suppression near defects^ - — 
In iron-based systems, where superconductivity exhibits 
both s-wave characteristics and a small coherence length, 
the situation is between the low-temperature and high- 
temperature superconductors. Consequently these sys- 
tems pose an intriguing question of how the effect of 
disorder on T c and the superfiuid density in these com- 
pounds compares to that in conventional BCS supercon- 



ductors and cupratesJ^ 

Recently, two irradiation experiments on Co-doped 
BaFe2As2 (Co-122) were performed to study the influ- 
ence of disorder . 11 ' 12 The temperature-dependent pene- 
tration depth measurements suggested an state, with 
strong nonmagnetic scattering in the unitary Hmitjei 
whereas transport measurements showed an s ++ state 
with weak scattering in the Born limits Both experi- 
ments showed a relatively small suppression of T c caused 
by nonmagnetic impurities induced by irradiation; these 
findings are consistent with an s ++ state, since super- 
conductivity with a sign changing order parameter is 
quite sensitive to nonmagnetic impurities »2ir— Reports 
in several iron-arsenide systems by different experimen- 
tal techniques are consistent with theoretical predictions 
of s wave, potentially nodal s-wave or sign reversing s- 
wave^ 

In this work we investigate the influence of random 
point defects introduced by proton irradiation on A(T) 
in Cao.5Nao.5Fe2As2 (CNFA) single crystals. We use the 
magnetic force microscopy (MFM) technique to deter- 
mine absolute values of A(T)^~— The CNFA single crys- 
tals, showing homogeneity, have been grown with a self- 
flux technique. Details of the sample preparation and 
characterization can be found elsewhere^ 



II. EXPERIMENT 



The 3 MeV protons are known to produce between one 
and a few tens of atomic displacements^ creating ran- 
dom point defects as well as nanoclusters with typical di- 
mensions of few nanometers. The CNFA sample was irra- 
diated with the total proton dose of 2xl0 16 cm -2 , which 
corresponds to an average distance (d) between defects 
of 2.8 nmpi The sample was cleaved, and its thickness 
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FIG. 1: (color online) Temperature dependent H C 2 along the c 
axis and within the ab plane in unirradiated Cao.5Nao.5Fe2 As2 
(black triangles), taken from Ref. 24, and proton irradiated 
Cao.5Nao.5Fe2As2 (red circles). 



measured to be around 28 yum, which is smaller than the 
penetration range of 40 /xm for the 3 MeV proton beam. 
Electrical resistivity in both unirradiated and irradiated 
samples were measured using a standard four-probe tech- 
nique. The sample was mounted in a rotatablc probe and 
measurements were performed in magnetic fields varying 
between and 9 T. MFM measurements described here 
were performed in a home-built low-temperature MFM 
apparatus^ Three samples, CNFA, irradiated CNFA 
(ICNFA), and a Nb reference film were loaded and in- 
vestigated in a comparative experiment within a single 
cool-down. The magnetic stray field calibration was per- 
formed by imaging vortices in a Nb reference as a function 
of applied magnetic fields Measurements of A were per- 
formed using the Meissner response technique »22r— The 
Meissner response curves were first measured as a func- 
tion of the tip-sample separation in the Nb reference with 
known A(T). Subsequently, the cantilever was moved to 
a sample of interest and the Meissner response curves 
were acquired. Direct comparison of measured curves 
yields the absolute value of A in a sample under investi- 
gation. Details of experimental technique are described 
elsewhere £2r— 

The reference Nb thin film (T c w 8.8 K) has a thickness 
of 300 nm and was grown by electron beam deposition. 
The T c of CNFA from transport measurements is 19.4 K 
and that of ICNFA is 17.8 K. The width of the super- 
conducting transition did not change after irradiation. 
No upturn in resistivity was observed at low tempera- 
tures, indicating that irradiation by protons results in the 
formation of nonmagnetic point-like scattering centers^ 
The MFM measurements were performed using a high- 
resolution Nanosensors cantilever— that was polarized 
along the tip axis in a 3 T magnetic field. Both Nb and 
CNFA samples were zero-field cooled for Meissner exper- 
iments; a magnetic field of a few Oe was applied above 
T c , followed by cooling for vortex imaging experiments. 



III. RESULTS 

A. H C 2(T) measurements 

Figure [T] shows the upper critical field H C 2(T) with 
H || c and H _L c (within the ab plane) for CNFA and 
ICNFA. H C 2(T) is linear in both samples and for both 

directions, with average slopes of ft ab = — d g^ 2 \ T = 
4 T/K and (3 C = -?§f-\ T = 2.2 T/K for CNFA, and 
pab = -^a^ = 3.8 T/K and 0° = -^| r<> = 2.3 
T/K for ICNFA. A modest superconducting anisotropy 

gab 

parameter 7 = = 1.85 — 1.65 for both CNFA and 
ICNFA samples points toward a three-dimensional be- 
havior. The superconducting coherence length £ can be 
expressed in the Ginzburg-Landau region as ^ql{T) « 
£0 / \/l — T/T c . In the case of a one- band model or two 
weakly coupled bands with similar Fermi surface prop- 
erties and pairing interactions the zero-temperature in- 
plane coherence length Q b and out-of-plane coherence 
length £q are given by the slope of the upper critical 
neldi^iS (£f) 2 w $ /2itT c (3 c and (£g) 2 « <i> /2nT c p ab . 
We obtain the zero-temperature Ginzburg-Landau val- 
ucs of £cW(0) - 2.8 nm and tf b CNFA (0) = 2.8 nm, 
which arc similar in magnitude to the short-coherence 
length cuprate and PuCoGa 5 superconductors. Within 
our measurement uncertainty no appreciable change of 
the coherence length took place after irradiation, al- 
though T c is suppressed by 10%. 

B. A(T) measurements 

Prior to measurements of the absolute values of A(T), 
vortex images were obtained under the same experimen- 
tal conditions for all samples. These measurements yield 
information about homogeneity of CNFA and ICNFA 
samples on a submicron scale (~ 100 nm). The well- 
formed single vortices in Nb and CNFA suggest the ho- 
mogeneity of the sample; however, the irregular shape of 
single vortex in ICNFA (elongated vortex in the diagonal 
direction of the image) suggests the presence of inhomo- 
geneity in the superfluid density on a sub-micron scale, 
which may be related to impurities introduced from irra- 
diation. We employed the following imaging procedure: 
First, a single vortex in the Nb sample was obtained at 4 
K after the stray field calibration of the MFM system^ 
Second, the MFM tip was moved on to CNFA and a 
single vortex image obtained, and third, a single vortex 
image was obtained after the tip was moved on to IC- 
NFA as shown in Figs. (Ha), (b), and (c). The line profile 
for each of the single vortices is shown in Fig. HJd). The 
intensity of the vortex center in different samples cor- 
relates with the magnitude of A, since all images were 
taken under the same conditions and with the same tip. 
Lower intensity corresponds to a larger A; therefore, A in 
ICNFA is much larger than that in the Nb reference. In 
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FIG. 2: (color online) Single vortex images in (a) the Nb ref- 
erence, (b) the unirradiated Cao.sNao.5Fe2As2, and (c) the 
irradiated Cao.sNao.5Fe2As2. (d) Comparison of single vor- 
tex profiles obtained from (a), (b), and (c). All images were 
obtained under the same experimental conditions in a single- 
cool down with the tip lift height of 300 nm at 4 K. The color 
scale bar refers to (a)-(c). 
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FIG. 3: (color online) Meissner response curves obtained from 
(a) the Nb reference (blue diamonds), (b) the unirradiated 
Cao.5Nao.5Fe2As2 (green diamonds), and (c) the irradiated 
Cao.5Nao.5Fe2As2 (red diamonds) at 4 K. The different slopes 
of the Meissner curves obtained from each sample indicate a 
systematic change of A. The inset: The Meissner curves for 
the unirradiated and irradiated samples are shifted along the 
horizontal axis to overlay the Meissner curve for the reference 
Nb sample. The difference of the penetration depths AA can 
be obtained from the values of the shift. 



addition, the magnitude of A among the superconducting 
samples can be inferred from the relative size of a single 
vortex: The larger the size, the larger is A. Therefore, A 
in ICNFA, showing the largest vortex size, is the biggest 
among them. 

To extract absolute values of A in ICNFA we performed 
the Meissner response measurements as described above. 
The Meissner curves as a function of the tip-sample sep- 
aration were obtained in all three samples, Nb, CNFA, 
and ICNFA (see Fig. [3J. The decay rate of the fre- 
quency shift 5f as a function of the tip-sample sepa- 
ration z provides the relative magnitude of A, i.e., the 
higher the rate d(Sf)/dz the larger the A. In bulk and 
thick films, the Meissner response force obeys a universal 
power-law dependence with tip-to-sample distance . 27 ' 28 
The force is given by Fueissner = A x /(A + z), where 
z is the tip-to-sample distance, A is a pre-factor con- 
taining information about the geometry of the magnetic 
tip, and f(z) ~ 1/z 3 . By shifting the f ICNFA (z) data 
with respect to distance in order to overlay it with the 
f Nb (z) curve, one can obtain the absolute values of 
Xicnfa{T) = \ N b(T)+A\{T), where AA(T) is the mag- 
nitude of the shift. The shift A A between the Nb and 
ICNFA data equals 320 nm, resulting in Xicnfa(0) = 
Ajvb(O) + AA(0) = 110 nm + 320 nm = 430 nm. Us- 
ing the same procedure we also obtained \cnfa{0) = 
X Nb (0) + AA(0) = 110 nm + 150 nm = 260 nm. Our 
experimental error is around 10% and depends on the 
magnitude of A and the system noise level. A key result 
of this work is that the A(0) values before and after ir- 
radiation differ significantly. This is in stark contrast to 
both the coherence length £, which shows little change 
after irradiation, as well as the small suppression in T c 
of 10%. The Meissner force MFM measurements of the 
ICNFA sample were performed after cleaving followed 
by irradiation. The ICNFA sample was remcasured after 
polishing. Both measurements showed the same A within 
experimental uncertainty This indicates that irradia- 
tion does not noticeably affect sample quality. There- 
fore we can neglect the degradation of the sample sur- 
face for Meissner screening currents. It should be noted 
that our parameter free method of using the Nb reference 
sample is based on the assumption of a universal scal- 
ing function F(z) for the Meissner force. This approach 
is valid for type-II superconductors, where the electro- 
magnetic response is local, i.e., n = A/£ 3> 1. Here we 
neglected higher order corrections in 1/k. Our Nb film 
has k = A/£ = 110nm/107iTO « 10. The large k value 
in Nb allows direct comparison of Meissner responses be- 
tween the Nb reference and CNFA, which results in good 
agreement by overlaying the Meissner curves, shown as 
insets in Fig. [3] and Figs. IDJa)-(b). Our novel method of 
using a reference sample is justified a posteriori because 
the Meissner curves would not overlay with one another 
just by shifting them. 

The temperature-dependent Meissner response curves 
measured in both CNFA and ICNFA samples are shown 
in Figs.HJa) and (b). The gradual variation of the Meiss- 
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FIG. 4: (color online) Temperature dependent Meissner 
curves for (a) the unirradiated Cao.sNao.5Fe2As2 and (b) the 
irradiated Cao.sNao.5Fe2As2 samples. Insets in (a) and (b) 
show overlaid temperature dependent Meissner curves at 4 
K, validating our procedure for extracting A(T). (c) Tem- 
perature dependent A(T) in both unirradiated and irradiated 
Cao.5Nao.5Fe2As2 samples determined from (a) and (b). (d) 
A(T) from (c) normalized by the T = value as a function of 
the normalized temperature. 



ner curves as a function of temperature indicates a sys- 
tematic change of A(T). The insets in (a) and (b) show 
the Meissner curves obtained at different temperatures 
but shifted to lie on top of the Meissner curve taken at 
T = 4 K; the curves overlay each other very well. The 
shift value for a given T to T = 4 K along the horizon- 
tal axis allows one to calculate A(T) at T. The result- 
ing X(T) and normalized A(T)/A(0) in both samples are 
shown in Figs. QJc) and (d), respectively. Results indi- 
cate that A(T) increases after proton irradiation; how- 
ever, the dependence of A(T)/A(0) on the normalized 
temperature T/T c is the same for both samples within 
our experimental uncertainty. The penetration depth ex- 
hibits the typical power-law behavior AA(T)/A(0) ~ T n 
with n 2 reported previously for doped iron-arsenide 
superconductors £ 



IV. DISCUSSION 

The radiation dose of 2xl0 16 cm -2 produced by 3 
MeV protons in a Cao.5Nao.5Fe2As2 sample causes the 
suppression of the superfluid density p s (0) « 1/A 2 (0) by 
about 60% whereas T c is only suppressed by 10%. We 
plot the value of the normalized p s (0) for ICNFA as a 
solid circle in the Uemura plot2& of disordered supercon- 
ductors in Fig. [SJ as well as theoretical results of one-band 
AG for d-wavc pairing (solid line) and two-band AG cal- 



culations for pairing (red open circles)^ Also shown 
are the BdG (Bogoliubov-de Gennes) calculations for d- 
wave pairing (red hatched circles), a Swiss cheese model 
far from the AG theory^ Our result bears similarity to 
the data for self-irradiated PuCoGag^S and He- irradiated 
YBCO high-tcmpcraturc superconductor, showing that 
T c is strongly immune to disorder relative to p s (0) ) 39 i 40 
contrary to the conventional AG theory for d-wave par- 
ing. By analogy we argue that the break-down of the 
AG theory is accounted for by the Swiss cheese model 
within the BdG lattice theory of short-coherence length 
superconductors^ which shows an abrupt suppression 
of the order parameter near point defects. This model 
describes the spatial dependence of the local density of 
states and the order parameter in the vicinity (within 
a few lattice constants) of a point-like nonmagnetic im- 
purity in the strong scattering limit, similar to holes in 
Swiss cheese. Franz and coworkers^ also reported the 
break down of the AG theory and strong suppression of 
p s (0) for d-wave paring. The effect is stronger in samples 
with small £/ao ratio (ao is the lattice constant). In the 
opposite limit, the AG theory is valid and the order pa- 
rameter is then suppressed uniformly in the entire sample 
because £ ^> ao,d, with d the average distance between 
impurities. In our sample the ratio of £o/ a o is approxi- 
mately 7, £ ~ 2.8 nm and d is about 2.8 nm, justifying 
the Swiss cheese scenario. 

It is worth noting that the T dependence of A(T) re- 
mains the same after irradiation as shown in Fig. 0] (d), 
while it changes in cuprates42r— This discrepancy may 
result from the nature of the multiband s-wave paring 
as well as the highly disordered nature of CNFA on the 
Ca/Na sites which lie above and below the iron layer. 
The fact that the temperature behavior of A(T) is ro- 
bust after irradiation may be ascribed to large intraband 
scattering with pairing and that the system itself is al- 
ready in the "dirty" limit prior to irradiation, consistent 
with its short coherence length and power-law depen- 
dence of A(T). Additional disorder (mostly in the iron 
layer) by proton irradiation therefore has little impact 
on the temperature behavior of A(T), while added inter- 
band scattering is detrimental to (increases) the absolute 
magnitude of A(0). 

The pair-breaking effect due to nonmagnetic scatter- 
ing in the AG theory can be quantitatively analyzed us- 
ing the normalized scattering rate in conjunction with 
A given by: g = hApo/ (2irkBT c op:o\Q), where Apo is 
residual resistivity change induced by irradiation, Apo = 
Po' r ~ Po mrr , T c q is the critical temperature before ir- 
radiation, and Ao is the penetration depth of the unir- 
radiated sampled The parameter g x and T c q are ex- 
pressed as ln(T c0 /T c ) = ^(1/2 + g x T c0 /(2T c )) - V(l/2), 
where ip(x) is the digamma function, based on the s 
scenario 4^ This pair-breaking result for T c is similar to 
that for conventional s-wave with magnetic impurities or 
d-wave with nonmagnetic impurities. Here the critical 
scattering rate parameter, where superconductivity van- 



ishes, is g 
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0.28 in the s pairing state. The 
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FIG. 5: (online color) Uemura plot of the superfluid den- 
sity in disordered short coherence length superconductors. 
T c o and p s o are values obtained from a pristine crystal; T c 
and p a are those measured after irradiation. The solid cir- 
cle represents the proton irradiated CNFA obtained in this 
work. For comparison we plot results of the one-band AG and 
BdG (Swiss cheese) calculations^ 7 - for d-wave pairing, two- 
band AG s ± -wave calculations 2 ^, and experimental results 
for self-irradiated PuCoGaj^ and helium irradiated YBCO 



V. CONCLUSION 



We reported the influence of random point disorder 
produced by proton irradiation on the superfluid density 
in Cao.5Nao.5Fe2As2. It leads to a dramatic change of 
A(0) after irradiation, in contrast to the small variation 
of T c and predictions by the AG theory. Both £(T) and 
A(T) show similar temperature behavior before and after 
irradiation. This behavior may be understood within the 
Swiss cheese model, the pair-breaking nature of s inter- 
band superconductivity, and a short coherence length, 
which considers the spatial dependence of the order pa- 
rameter and its strong suppression near defects at the 
atomic scale. Finally, the extracted normalized scatter- 
ing rate, in conjunction with the absolute value of A(T), 
is much larger than the critical scattering rate for the 
s ± pairing, confirming the break-down of the AG theory 
in these disordered superconductors. Further detailed 
multiband BdG model calculations combined with sys- 
tematic doping and irradiation studies may shed light on 
the suppression of superconductivity in this large class of 
iron-based superconductors. 



extrapolated critical scattering parameter, obtained us- 
ing Ap =30 p,Q cm and A =260 nm, is « 3.7. This 
value is much larger than that expected in the sce- 
nario, quantifying the break-down of the AG theory in 
irradiated iron-arsenide superconductors, where the ap- 
proximation of the uniformly impurity-averaged Green's 
function is not valid. Similar results were reported in 
Ba(Fei_ a: Co a ;)2As2 irradiated by protons^ and illustrate 
the generality of the Swiss cheese model for pair-breaking 
in this large class of high-temperature superconductors. 
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